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\Q ' Many types of stars show periodic variations of radius and brightness, which 

are commonly referred to as 'stellar pulsations'. Observed pulsational character- 
istics are determined by fundamental stellar parameters. Consequently, investi- 

Q\ • gations of stellar pulsations provide a unique opportunity to verify and refine our 

understanding of the evolution and internal structure of stars. However, a key 

OS ■ boundary condition for this analysis - precise information about the geometry of 

pulsations in the outer stellar envelopes - has been notoriously difficult to secure. 
Here we demonstrate that it is possible to solve this problem by constructing an 
■ 'image' of the pulsation velocity field from time series observations of stellar spec- 

tra. This technique is applied to study the geometry of nonradial pulsations in a 
i-, ■ prototype magnetic oscillating (roAp) star HR3831. Our velocity map directly 

demonstrates an alignment of pulsations with the axis of the global magnetic 
field and reveals a significant magnetically induced distortion of pulsations. This 

J^j ■ observation constitutes a long-sought solution of the problem of the pulsation 

geometry of roAp stars and enables very stringent tests of the recent theories of 
stellar magneto-acoustic oscillations. 

Subject headings: line: profiles — stars: chemically peculiar — stars: imaging — 
stars: individual (HR3831) — stars: oscillations 



1. Introduction 

A classical solution (Unno et al. 1989) of hydrodynamic equations expresses monope- 
riodic nonradial pulsational fluctuations of an idealized spherically symmetric star in terms 
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of a spherical harmonic function. Approximations of this simple stellar pulsation theory are 
adequate for a significant fraction of oscillating stars and, therefore, historically investiga- 
tions of stellar pulsation geometry were limited to assigning spherical harmonic 'quantum 
numbers' i and m to each observed pulsation mode. This problem is solved by several 
varieties of the mode identification technique (Aerts, De Pauw, & Waelkens 1992; Watson 
1988), which, with different degrees of success, are applied to photometric and spectroscopic 
observations of pulsating stars. However, pulsation modes in real stars are often significantly 
distorted by rapid stellar rotation (Lee & Saio 1990) and global magnetic field (Dziembowski 
& Goode 1985; Bigot & Dziembowski 2002; Saio & Gautschy 2004). These effects result in a 
considerably more complicated pulsation geometry, irreducible to a single spherical harmonic 
function. Up to now no objective technique was available for a comprehensive and reliable 
analysis of such distorted nonradial pulsations, but it was understood that a detailed study 
of high-resolution stellar spectra is the most promising approach to resolve these difficulties. 

In general, all inhomogeneities in the outer layers of stellar envelopes, including pul- 
sational disturbances, produce distortions in the stellar line profiles. As a star rotates 
and starspots appear and disappear from the visible hemisphere, spectral line shapes and 
strengths are modulated periodically. Doppler imaging (DI, e.g., Vogt, Penrod, & Hatzes 
1987) is a powerful technique to decipher rich information about the physics and geometry 
of the surface structures contained in the stellar spectrum variability. The Doppler mapping 
method is based on the resolution of the stellar surface provided by the effect of stellar rota- 
tion on the shapes of absorption spectral lines. If the projected rotational velocity of a star 
is large enough for the shapes of spectral lines to be dominated by the rotational Doppler 
broadening, there exists a well-defined mapping relation between distortions in a line profile 
and positions of the corresponding inhomogeneities on the stellar surface. Thus, a single 
Doppler-broadened line profile represents a one-dimensional projection of the visible stellar 
hemisphere, with all the structure blurred along the lines of constant Doppler shift. Spec- 
troscopic observations at different rotation phases can be used to monitor the evolution of 
dips and bumps in line profiles, thus obtaining a set of non-redundant projections of a stellar 
surface map, and eventually to reconstruct a two-dimensional image from a time series of 
line profile observations. This indirect stellar surface imaging method was first applied to 
the problem of mapping chemical inhomogeneities in early-type stars (Khokhlova, Rice, & 
Wehlau 1986) and was subsequently extended to investigation of temperature spots (Vogt et 
al. 1987) and complex magnetic topologies (Semel 1989) in active late-type stars and global 
magnetic fields (Piskunov & Kochukhov 2002) in chemically peculiar stars. Vogt & Penrod 
(1983) have emphasized that rotational broadening of stellar spectra also aids in probing 
the horizontal structure of nonradial pulsations. However, subsequent spectroscopic studies 
of rapidly rotating oscillating stars (Gies & Kullavanijaya 1988; Kennelly, Walker, & Mer- 
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ryfield 1992; Hao 1998) used only resolution of pulsation structures along stellar longitudes 
and restrictively approximated velocity field with a single spherical harmonic. 

In a recent development Berdyugina, Telting, & Korhonen (2003) presented a new ex- 
ploration of the capabilities of the DI technique in its application to mapping of stellar 
nonradial pulsations. These authors employed a temperature DI code in modeling pulsa- 
tional line profile variations and assumed that pulsational fluctuations at the stellar surface 
can be represented by a rotation of a fixed pattern. As emphasized by Kochukhov (2004), 
the latter assumption restricts pulsational mapping to the reconstruction of a spherical har- 
monic geometry, which is aligned with the stellar rotation axis and is characterized by a 
unique non-zero azimuthal number m. This condition is not fulfilled for zonal modes and 
is invalid for all types of oblique stellar pulsations, hence the method of Berdyugina et al. 
(2003) cannot be applied to study oscillations of roAp stars. 

Here we present the first application of an alternative technique. In our approach 
the principle of Doppler imaging is extended to the reconstruction of the time-dependent 
velocity field. The temporal and surface variation of the velocity vector are represented with 
a superposition of the two constant surface distributions: 

V(t, 9, 0) = V c (6, 0) cos(wt) + V s (6, 0) sm(ut), 

where u is the pulsation frequency and 9, <fi are usual spherical coordinates on the stellar 
surface. The V c and V s vector maps are recovered directly from the observed line profile 
variability without imposing any specific global constraints on the pulsation geometry. This 
is equivalent to mapping a two-dimensional surface distribution of the pulsation amplitude 
and phase for each velocity component. 

The foundations of the pulsational Doppler mapping and description of its computer 
implementation were presented by Kochukhov (2004). We refer the reader to this paper for a 
detailed explanation of the technique and discussion of the numerical simulations which were 
used to evaluate performance and intrinsic limitations of the new surface mapping method. 

2. Pulsational Doppler mapping of HR 3831 

The pulsation Doppler inversion is applied to time-resolved observations of the Nd III 
A 6145 A line in the spectrum of the well-known roAp star HR3831 (HD 83368). A total of 
1860 spectra of this object were obtained over the period of 11 nights using the Coude Echelle 
Spectrograph fiber-linked to the Cassegrain focus of the 3.6-m telescope at the European 
Southern Observatory. A preliminary discovery report of the pulsational variations in indi- 
vidual spectral lines in HR3831 using these observational data was published by Kochukhov 
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& Ryabchikova (2001). Our observations evenly sample the 2 < ?851976 rotation period of the 
star (Kurtz et al. 1997) and are characterized by the resolving power of A/AA = 123,000 
and the signal-to-noise ratio of 110-160 pixel^ 1 . A 70 s exposure time was chosen to ensure 
an appropriate sampling of the 11™67 oscillation period. Fig. 1 shows an example of the 
rapid line profile variation of HR3831 and presents a comparison between the observed and 
calculated spectra. 

Reconstruction of the pulsation velocity field of HR3831 took into account an inho- 
mogeneous surface distribution of neodymium which significantly distorts the mean Nd III 
line shapes. The maps of chemical abundance and pulsational fluctuations were recovered 
simultaneously in a self-consistent manner. The resulting Nd abundance map agrees very 
well with the distribution derived by Kochukhov et al. (2004) using a different DI code. 

Nonradial pulsations in roAp stars are characterized by very small photometric ampli- 
tudes. The respective pulsational temperature variation does not exceed ~ 10 K (Medupe 
& Kurtz 1998). Such temperature changes do not have a noticeable effect on the stellar line 
profiles and can be safely neglected in our modeling. 

The formation region of the Nd III A 6145 A line is located high in the atmospheres 
of roAp stars (Ryabchikova et al. 2002). The horizontal motion due to the high-overtone 
p-mode pulsation is expected to be either negligible (Bigot & Dziembowski 2002) or consid- 
erably smaller (Saio & Gautschy 2004) than the radial motion at this atmospheric height. 
Trial mapping of the full vector pulsation velocity distribution demonstrated that, whereas 
including the horizontal pulsation motion substantially increases the number of free param- 
eters, it does not lead to an appreciable improvement of the fit to observations. Taking this 
into account we limited the pulsation DI reconstruction to the vertical velocity component 
which is represented by the Vf and velocity amplitudes. 

The Doppler maps of the pulsation velocity field of HR3831 revealed with our technique 
are illustrated in Fig. 2. The Vf map shows a clear oblique dipolar pattern with the local 
pulsation amplitude reaching up to 4.0 kms -1 . On the other hand, the amplitude of the 
stuctures in the second map, V r s , does not exceed a few hundred ms" 1 . Such a relation of 
the two maps is typical of nearly axisymmetric pulsations (see Sect 5.5 of Kochukhov 2004). 

3. Discussion 

Indirect imaging of the surface structure of nonradial oscillations has important im- 
plications for stellar astrophysics. Results of the first Doppler inversion of the pulsation 
velocity presented here can be regarded as an ultimate solution of the long-standing problem 
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of the pulsation geometry of roAp stars. Until now the high-overtone p-mode oscillations 
in these unique objects were investigated primarily with the high-speed photometric tech- 
niques. Analyses of light variations of roAp stars (Kurtz 1990) and, in particular, HR3831 
(Kurtz et al. 1997), demonstrated that the amplitude and phase of pulsations are modulated 
with the stellar rotation. A coincidence of the times of maxima of pulsation amplitude and 
magnetic field extrema was interpreted within the oblique pulsator framework (Kurtz 1982). 
This phenomenological model postulated that pulsations are represented by an axisymmetric 
dipolar (£ = 1) mode inclined with respect to the stellar rotation axis and aligned with the 
axis of the global magnetic field. Subsequent theoretical studies were carried out by Shiba- 
hashi & Takata (1993), Takata & Shibahashi (1995) and, most recently, by Saio & Gautschy 
(2004) assuming a dominant role of the magnetic field. These calculations suggested that 
important deviations of the mode geometry from the purely dipolar topology might occur, 
although pulsation structure was still expected to remain almost axisymmetric and aligned 
with the magnetic field. A totally different picture has emerged from the theoretical inves- 
tigation by Bigot & Dziembowski (2002) who focused on the essential role of the centrifugal 
force resulting from the stellar rotation. According to their theory pulsations in roAp stars 
are described by a superposition of all £ = 1 harmonic components and are neither axisym- 
metric, nor necessarily aligned with the magnetic field. Limited information content of the 
photometric observations did not allow astronomers to distinguish between these competing 
approaches to the roAp pulsation structure. 

The surface velocity mapping technique is not limited to the spherical harmonic param- 
eterization of pulsational fluctuations and can be used to explore the physics of nonradial 
oscillations in rotating magnetic stars. In particular, it is able to supply a crucial obser- 
vational constraint required for testing the alternative theories of roAp pulsations. The 
pulsation map (Fig. 2) reconstructed here for the prototype roAp star HR3831 displays an 
obvious correlation with the stellar magnetic field geometry. Kochukhov et al. (2004) found 
that the dominant dipolar component of the global field in HR3831 is inclined by about 
90° relative to the stellar rotation axis and is oriented at rs 180° in longitude. At the same 
time, the V r c velocity map presented in Fig. 2 indicates that pulsation geometry of HR3831 
is dominated by an axisymmetric structure. Furthermore, it is clear that the two maxima 
of the pulsation amplitude coincide with the magnetic poles. This is the first independent 
verification of the alignment of pulsations and magnetic field in oscillating magnetized stars. 

Our pulsational mapping of HR3831 suggests a dominant role of the magnetic perturba- 
tion of the p-modes and considerably less important influence of the stellar rotation. Indeed, 
the pulsation Doppler image of the surface velocity field reveals no significant asymmetry 
which can be attributed to the oblique £ — 1, \m\ = 1 components. Therefore, we conclude 
that pulsations in HR3831 cannot be described by the theory brought forward by Bigot 
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& Dziembowski (2002), who attempted to fully take into account the rotationally induced 
mode perturbations but limited their calculations to a relatively weak magnetic field. 

Fig. 3 illustrates an analysis of the dominant axisymmetric pulsation structure recon- 
structed with our technique. This figure shows the latitudinal dependence of the Vf pulsation 
amplitude in the reference frame of the stellar magnetic field. The inferred latitudinal trend 
is clearly different from the one expected for a pure oblique i = 1 pulsation. Axisymmetric 
multipolar decomposition of the observed latitudinal dependence suggests that oscillations 
in HR3831 can be roughly described with a superposition of t = 1 and t = 3 harmonic 
components. The octupolar component has the same sign as the dipolar contribution and 
an amplitude smaller by a factor of two (A^ =3 /A^ =1 =0.50±0.04). In other words, we discover 
that pulsations in HR3831 are strongly confined to the magnetic field axis. This pulsation 
geometry is in remarkable agreement with the predictions obtained in recent calculations as- 
suming the dominant role of the magnetic field (Saio & Gautschy 2004). This concordance of 
the empirical model and theory documents a significant breakthrough in our understanding 
of the physics of p-mode oscillations in magnetic stars and possibly confirms an important 
role of damping by the magnetic slow waves investigated by Saio & Gautschy (2004) in the 
context of excitation of the roAp-type pulsations. 

In addition to revealing the true pulsation geometry of roAp stars, our precise velocity 
map recovered for HR3831 opens numerous possibilities for asteroseismological determina- 
tion of various stellar parameters. For instance, Takata & Shibahashi (1995) argued that 
the presence of a large quadrupolar component in the magnetic field geometry is expected to 
induce £ = and £ = 2 axisymmetric structures in the pulsation velocity distribution. How- 
ever, Fig. 3 convincingly demonstrates that the even multipolar components are very weak 
in HR3831. Hence, the magnetic field of this star must be basically dipolar - a conclusion 
obtained independently from an analysis of any classical magnetic field observables. 

This work was supported by the Lise Meitner fellowship granted by the Austrian Sci- 
ence Fund (FWF, project No. M757-N02). We acknowledge extensive use of the SIMBAD 
database, operated at CDS, Strasbourg, France. 
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Fig. I. — Observed variation of the Nd III A 6145 A line (histogram) is compared with 
the best fit synthetic spectra (solid line). The left panel shows the observed and synthetic 
spectra for 20 phases (out of 1860 used in the inversion) close to the maximum of pulsational 
fluctuations (rotation phases tp = 0.966-0.973 according to the ephemeris of Kurtz et al. 
1997). The right panel presents the difference between individual profiles and the average 
line shape. The spectra for consecutive pulsation phases are shifted in the vertical direction. 
The bars at the bottom of the panels show the vertical scale in units of the continuum flux. 
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Fig. 2. — Results of the first Doppler imaging reconstruction of the stellar pulsation velocity 
field. Rectangular panels show the Vf (left) and V* (right) vertical velocity amplitude maps. 
The contours of equal velocity are plotted over the greyscale images with a step of 1.0 
kms -1 . The horizontal dotted line shows the lowest visible latitude corresponding to the 
angle between the stellar axis of rotation and the observer's line of sight ('inclination angle', 
i = 68°) adopted in the analysis of HR3831. The plots below present spherical projection 
of the respective pulsation velocity maps. The star is shown at 4 different aspect angles 
corresponding to the rotation phases tp = 0.00, 0.25, 0.50, and 0.75. The grid at the stellar 
surface is plotted with a 30° step in longitude and latitude. 
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Fig. 3. — Analysis of the latitudinal dependence of the V r c pulsation velocity amplitude. The 
symbols show average velocity amplitude as a function of the cosine of the angle 9 counted 
from one of the pulsation poles. The solid line illustrates the fit with a superposition of 
axisymmetric £ = 0, 1, 2, and 3 spherical harmonic components. Contributions of individual 
multipolar components are shown by the dashed lines for £ = 1-3 and by the dotted line for 
£ = 0. 



